UNCLASSIFIED 


AD  NUMBER 


AD492301 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release^  distribution 
unlimited 


EROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  MAY  1950. 
Other  requests  shall  be  referred  to  US 
Navy  Office  of  Naval  Research^  Washington^ 
DC. 


AUTHORITY 


ONR  notice,  27  Jul  1971 


THIS  PAGE  IS  UNCLASSIEIED 


NOTICE;  l^en  grtvwnwiii'it  or  other  dmrlags^  ipeei- 
fleetlodf  or  other  date  ere  used  for  say  purpose 
other  than  In  conneetloa  vlth  e  definitely  nleted 
gorenaunt  proeuriunt  operetion,  the  U.  8. 
Oorenaent  thereby  Incurs  no  responsibility^  nor  say 
obligation  uhetsoeverj  sad  the  feet  that  the  Qovem- 
sant  nay  have  fozaulated,  fUrnlshe^  or  In  ai^^  way 
supplied  the  said  drawings^  speelfleatlcns,  or  other 
data  Is  not  to  be  regarded  by  lapUeatlon  or  other- 
vise  as  In  any  Manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rl^ts 
or  pemloslon  to  aaauflaeture,  use  or  sell  ai^ 
patented  Invention  that  nay  la  any  vay  be  related 
thereto. 


Best 

Available 

Copy 


UNIVERSITY  OF  CALIFORNIA 


DEPARTMENT  OF  ENGINEERING 
LOS  ANGELES 


INVESTIGATION  OP 

Atmospheric  Diffusion  Processes 
BY  means  op 

Experimental,  analytical,  and  numerical 


J.  ^OPPENDIEK 

1.  t.  Vehrencamp 


Work  Sponsored  by  Office  of  Naval  Research 
(Geophysical  Branch) 


L.  M.  K.  BOELTEH.  QiAIRMAN  OF  THE  DEPARTMENT  OF  DfCINBERIMO 


i 


TABLE  OF  CONTENTS  . 


Suomry 

Nonanclatiire 

Introdttotlon 

Basic  Ldd3'  OirfUslon  Relations 

Experimental  Techniques  of  Determining  Eddy  Diffusivities^ 
Thermal  Conductances,  and  Drag  Coefficients 

A,  Convective  Heat  Tlransfer 

B.  Momentum  Transfer 

Analytical  cvA  Numerical  Eddy  Diffusion  Analyses 


A.  Analytical  Analyses 

B. .  Analyses 


A*  Eiqperlmental  Determination  of  Convective  Heat  Plows, 

Convective  Ccmductances,  Shear  Stresses,  Drag  Coefficients, 
and  Eddy  Oiffusivities 

1*  Heat  Transfer 

2,  Komentun  Transfer 

B*  Numerical  Eddy  Diffusion  Analyses  of  Transient  Air 
Temperature  Data 

C.  Eddy  Diffusion  Analyses  Found  in  the  Literature 
DISCUSSION 

ACKNOWLBDOIIENTS 

APPENDIX 

\ -  -  ■'  ^ 

REFERENCBS  _ 


1 

2 

6 

6 

10 

10 

12 

16 

16 

M _ 

2S 

2d 

2d 

29 

37 

41 

45 

49 

50 

52^  - 


!  i 

I  1 

i 

T 

I 

I 

'  I 

i 

I 

j 


i 


4 


I 


1-/ 

-■nui  - 


paper  concerns 


SUmRY 


an  investigation  of  the  atmospheric  diffusion 


processes  by  means  of  experimental,  analytical,  and  numerical  methods. 

Experimental  momentum  transfer  and  heat  transfer  analyses  are  presented. 

Shear  stresses  in  the  surface  layer  are  measured  directly  by  means  of  a  shear 
meter;  momentum  eddy  diffusivities  and  drag  coefficients  over  flat  ground  are 
determined  from  shear  stress  and  wind  velocity  data.  Convective  heat  transfer 
rates  in  the  surface  layer  are  measured  indirectly  by  means  of  a  heat  flow  meter 
and  a  radiometer;  ccnvectlve  conductances  are  determined  from  the  convective  heat 
flow  and  air  temperature  data. 

Two  new,  analytical,  periodic,  convective  heat  flow  solutions  for  the  atmos¬ 
pheric  system  are  derived.  One  solution  pertains  to  an  eddy  diffusion  system 
in  which  the  boundary  temperature  varies  sinusoidally^  with  time  and  the  eddy  dlf- 
fusivity  varies  sinusoidally  with  time  but  is  independent  of  height.  The  other 
solution  pertains  to  an  eddy  diffusion  system  in  which  the  boundary  temperature 
varies  sinusoidally  with  time  and  the  eddy  diffusivity  varies  sinusoidally  with 
time  and  linearly  with  height. 

A  simple  heat-momentum  transfer  analogy  for  the  surface  layer ^ws  been  devel¬ 
oped  which  relates  some  of  the  pertinent  heat  and  momentum  transfer  variables  in 
the  atmospheric  diffusion  s; stem.  V 

Three  numerical  methods  of  deteimnlng  eddy  diffusivities  as  a  function  of 
time  and  height  from  transient  air  tmmperature  and  humidity  profiles  are  presented. 
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coefficients  of  the  Fourier  cosine  series,  equation  (22),  ®P 

2 

shear  meter  surface  area  exposed  to  air  flow,  ft 

A 

constant  in  equation  (31) »  ft  Ar 
constr.nt  in  equation  (31),  ftA** 
constant  in  equation  (31),  dimensionless 
function  of  ^  defined  in  ecuation  (42) 
constant  in  equation  (24) 
coefficients  in  eouation  (25),  ®P 
water  vapor  concentration, 
drag ■ coefficient 
constants  in  equation  (14),  ft^A** 
fluid  heat  capacity,  Btu/i*  °P 


constant  in  ecuation  (47) 
coefficients  in  eouation  (43),  °P 
base  of  natural  logarithms 

unit  thermal  convective  conductance,  BtuAr  ft^ 
drag  force  on  the  shear  meter,  # 
acceleration  of  gravity,  ftA** 

coefficients  in  the  Fourier  cosine  series,  ecuation  (39), '**F 
Hankel  function  of  the  zero  order 

nTT 

ratio  of  the  vertical  eddy  diffusivity  for  mass  transfer  to 
the  vertical  eddy  diffusivity  fo*.*  heat  transfer 

-  thersial  conductivity  of  air,  BtuAr  - r 

spring  constant,  Ibs/inoh 

Karman  constant 

latmt  heat  of  >raporiiation  of  water,  Btu/# 
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modulus  of  the  complsx  function 
positive  integer 
constant  in  equation  (43) 
constants  In  equation  (4d)« 

vertical  convective  heat  transfer  rate  per  unit  area,  Btu^r  ft^ 

vertical  convective  heat  transfer  rate  per  unit  area  at  the 
ground,  Btu/hr  ft^ 

heat  transfer  rate  per  unit  area  due  to  evaporation  or  condm- 
satlon  of  vater  vapor  at  the  interface,  Btu^u*  ft^ 

heat  transfer  rata  per  unit  area  due  to  radiation  from  water 
vapor,  CO2,  and  dust  in  the  atmosphere,  BtuA** 

heat  transfer  rate  per  unit  area  into  or  out  of  the  ground  at 
the  Interface,  Btu/nr  ft^ 

groxuid  radloslty  per  unit  area  (emitted  and  reflected  radiation), 
BtuAr  ft< 

solar  irradiation  per  unit  area,  Btu/hr  ft^ 

total  hemispherical  radiation  per  unit  area  at  the  ground, 

Btu^  ft2 

a  variable  defined  by  equation  (40),  ft/(hr)i 
a  function  of  r  defined  by  equation  (42) 
an  angle  whose  tangent  is  described  for  equation  (50) 
potential  temperature,  equal  to  T  Pi,  ^ 
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amplitude  of  the  sinusoidal  boundary  tenqperature  wave,  ^ 

mean  air  temperature  at  the  laminar  sublayer  >  turbulent 
layer  Interface  (see  Figure  15),  ^F 

mean  air  temperature  at  the  reference  height,  v*  ^ 

mean  air-earth  interface  temperature, 

mean  fluid  velocities  in  the  x,  y,  and  s  directions,  .respec¬ 
tive:^,  ft^ 

mean  air  velocity  at  the  laminar  sublayer  -  turbulent 
layer  interface,  ftAr _  _ 

mean  air  velocity  at  the  reference  height,  ,  ftAr 


X.  y,  X  the  oarteeian  coordinates  (the  earth's  surface  is  in  the  Jff 
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‘v  V  S 


,  V  *xi  S 


laminar  sublayer  thickness^  ft 
rouj^hness  parameter,  ft 
a  function  of  z  .defined  by  equation  (23) 


«;rey  body  earth  surface  emlsslvity 

ratio  of  the  vertical  eddy  dlffusivity  for  heat  transfer  to  the 
vertical  eddy  dlffusivity  for  momentum  transfer 

a  variable  defined  by  equation  (33}«  hr 

adiabatic  lapse  rate,  ^P/ft 

air  density, 

spring  deflection,  inches 

thermal  eddy  dlffusivitles  in  the  x,  y,  and  z  direetionsj 
ft^/hr 

vertical  eddy  dlffusivity  for  mass  transfer,  ft^A** 

vertical  eddy  dlffusivity  for  momentum  transfer,  ft^Ar - - 

coefficients  in  eouatlons  (52),  (55),  and  (56) 
time,  hr St 

period  of  the  sinusoidal  boundary  temperature  and  dlffusivity 
vaves,  hrst 

A 

a  variable  defined  by  equation  (16),  ft 
a  variable  of  Integration  in  equation  (22) 
a  function  of  X  defined  in  equation  (23) 
absolute  air  viscosity,  hr/ft^ 
constant  in  equation  (24),  l/ft^ 
coefficients  in  equation  (25),  l/tt^ 
klnmatic  viscosity,  ft^Ai* 
mass  density  of  fluid,  # 

Stefan-Boltsmann  constant,  17t3X10"^®  Btu/ft^(®H)^ 
fluid  ^rtMsr  stress,  ~ 

fluid  shear  stress  at  the  ground,  #/ft^ 
amplitude  of  the  complex  function 


waves,  ^undai*y  tenparature  and  diffusivlty 


Dimensionless  ffadnH 


IMTRODUCTIOH 


A  iiuniber  of  transfer  processes  -which  are  associated  with  current  atmos¬ 
pheric  diffusion  problems  are  the  diffusion  of 

1)  air  pollutants  from  indust'-i.al  and  other  sources « 

2)  water  vapor  from  lakes « 

3)  water  vapor  from  snow  banks* 

U)  heat  and  momentum  from  orchards*  and 

5)  smokes  and  poisonous  gases  from  grenades  or  bon^s. 

Each  of  these  transfer  processes  falls  into  the  general  category  of  atmospheric 
diffusion. 


The  diffusion  of  heat*  mass*  and  moiaentum  in  the  atmosphere  is  achieved 
by  a  complicated  turbulence  mechanism  which  at  present  does  not  appear  to  be 
thoroughly  understood.  An  exact  fxmdamental  relation  describing  eddy  diffusion 
has  not  yet  been  determined.  Several  approximate  eddy  diffusion  relations  are 
available*  however.  One  of  the  most  useful  of  these  approximate  diffusion 
relations  is  the  one  which  relates  the  rate  of  transfer  of  a  quantity  to  a 
potential  gradient  and  an  eddy  dlffusivity. 

The  eddy  diffusion  rate  of  heat,  mass*  or  momentum  at  the  earth's  surface* 
in  particular,  may  be  expressed  in  terms  of  flow  potential  differences  and 
corresponding  transfer  conductances  and  coefficients.  For  example*  vertical 
convective  heat  transfer  rates  at  the  earth's  surface  can  be  e^qpressed  in 
terms  of  air  temperature  differences  and  unit  thermal  convective  conductances* 
and  vertical  momentum  transfer  rates  at  the  earth's  surface  oan  be  expressed 
in  terms  of  wind  velocities  and  drag  coefficients. 

If  quantitative  atmospheric  turbulence  data  such  as  eddy  dlffusivity 
profiles,  unit  convective  conductances*  or  drag  coefficients  are  available, 
it  is  possible  to  estimate  heat*  mass*  and  monentua  transfer  rates  or  flow 
potcnitlal  concentrations  in  the  atmosphere  for  point  source*  lino  soiuree*  or 
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area  source  diffusion  systems.  These  basic  turbulent  diffusion  data  must  be 
available  before  satisfactory  solutions  to  some  of  the  current  atmospheric 
diffusion  problems  can  be  effected. 

The  follovrlng  paragraphs  consist  of  a  discussion  of 

1)  some  of  the  basic  eddy  diffusion  relatims  to  be  considered  t 

2)  some  experimental  methods  of  determining  eddy  dlffusivities,  thermal 
conductances,  and  drag  coefficients, 

3)  analytical  and  nunerlcal  eddy  diffusion  analyses,  and 

It)  the  analysis  of  several  sets  of  mlcrometeorologleal  measurements  by 
the  methods  outlined  in  this  paper. 
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BASIC  EDDY  DIFFJSIOM  DELATIONS 

The  eddy  transrer  rate  of  a  quantity  in  the  absence  of  atnospherie  ther¬ 
mals  is  of-en  expressed  in  teruir  of  an  eddy  diffuaivity  and  a  potential  gradient. 
For  example*  the  vertical  rate  of  convective  heat  transfer  is*  (Reference  1) 


(a). 


•.(It  » '■) 


(1) 


where* 


y  •  air  density 

■  air  heat  capacity 
•  vertical  thermal  eddy  diffusivity 
T  '  mean  air  temperature 
r  ■  adiabatic  lapse  rate 


Similarly*  the  vertical  rate  of  momentum  transfer  in  the  absence  of  atmospheric 


1^' 

k 

I 


f 
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thermals  expressed  in  tei'ius  of  the  fluid  shear  stress  is*  (Reference  2) 

r  -  p€  /|Ji\  (2) 

•jto«  \3x/ 

iriiere* 

p  >  mass  density  of  air 

«_  •  moment\ira  eddy  diffusivity  in  the  a  direction 

*lfo» 

U  mean  air  velocity 

A  mass  transfer  equation  similar  to  equations  (1)  and  (2)  can  also  be  written. 

If  a  heat  balance  and  a  force  balance  are  made  on  a  differential  lattice 
in  a  twbulent  flow  system*  the  transient  eddy  diffusion  equations  for  heat 
arid  momentum  transfer  can  be  derived*  respectively.  For  example*  the  heat 
transfer  diffusion  equation*  when  the  molecular  conduction  terms  are- small* 
is*  (Reference  3) 


> 


y 


iiheret 

6  •  time 

I 

U.  V.  «  mean  fluid  velocities  In  x,  y,  and  z  dlrectlonst  respectively  .  , 

c  .  «  ,  e  ■  thermal  eddy  diffusivities  in  x,  y,  and  s  directions^  I 

1^11  respectively  J 

The  Navler>Stokes  equations  irhlch  were  modified  by  Osborne  Reiynolds  to  include  | 

! 

the  fluctuating  velocity  components  are  similar  to  the  heat  transfer  diffusion  j 

equation  given  above  except  for  the  presence  of  a  pressure  gradient  term  in  each 
of  the. three  hydrodynamic  equations  (Reference  U)* 

Sometimes  it  is  convenient  to  express  vertical  heat  and  momentum  transfer 
rates  at  the  earth’s  surface  (the  boundary  equations)  in  terms  of  flow  potential 
differences  and  corresponding  transfer  conductances  and  coefficients  rather  than 
flow  potential  gradients  and  eddy  diffusivities.  This  procedure  has  long  been 
used  in  the  field  of  engineering  in  connection  with  fluid  flow  and  heat  transfer 

systems  of  finite  dimensions.  The  vertical  rate  of  convective  heat  transfer  at _ [ 

the  earth's  s'Jrface  can  be  expressed  as  .  •  ^ 


(4) 


where, 

■  unit  thermal  convective  conductance 
»  mean  air-earth  interface  temperature 
>  mean  air  temperature  at  a  reference  height 
The  shear  stress  can  be  expressed  as 

\  ;  (5) 

where, 

«  drag  coefficient 

mean  air  velocity  at  a  reference  height 
The  basic  eddy  diffusion  relations  that  have  been  noted  above  are  required 
for  the  diffusion  analyses  which  are  presented  in  the  remainder  of  this  paper. 
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MPSRRtEMTAL  TECHNIQUES  OF  DBO^ERMININQ 
EDDY  DIFFUSIVITI5S,  THERMAL  COrmUCTAHCBS,  AND  DRAG  COEFFICIENTS 

A.  Vertical  Convective  Heat  Transfer. 

Because  of  the  recent  development  of  two  simple  thermal  instrviments,  it 

is  possible  to  measure  vertical  eddy  diffuslvity  profiles  and  unit  thermal 

conductances.  A  heat  flow  meter  (References  5  and  6)  has  been  developed  at 

0 

the  University  of  California  which  can  be  used  to  measure  heat  flow  rates  into 
or  out  of  the  earth's  surface.  This  heat  flow  meter  consists  of  laminated 
sheets  of  thin  bakelite  with  an  embedded  thermopile.  Heat  flow  through  this 
meter  is  measured  in  terms  of  the  voltage  response  of  the  calibrated  thermo¬ 
pile.  A  total  hemispherical  radiometer  (Reference  7)«  also  developed  at  the 
University  of  California,  can  be  employed  to  measure  the  total  solar  and 
nocturnal  irradiation  upon  the  earth's  surface  as  well  as  the  radiosity  of 
the  earth's  surface.  Briefly,  this  instrument  consists  of  a  horiaontal  heat 
meter  idiose  upper  surface  is  blackened  and  whose  lower  surface  is  surfaced 
with  a  sheet  of  aluminum.  Both  the  upper  and  lower  surfaces  of  the  heat 
meter  are  exposed  to  air  streams  of  equal  velocity  originating  from  a  small 
blower.  If  a  complete  heat  rate  balance  is  made  on  this  aystem,  it  'can  be 
shoim  that  the  total  hemispherical  radiation  falling  upon  the  horizontal  heat 
meter  surface  is  equal  to  a  constant  times  the  voltage  drop  across  the  heat 
meter  thermopile  plus  a  datum  term. 

If  a  heat  balance  is  made  at  the  earth's  surface,  the  convective  plus 
evaporative  heat  transfei*  rates  may  be  expressed  as^ 


where,  _ 

a  convective  heat  transfer  rate  per  unit  area  at  the  ground 

A I  oofiv^ 

The  signs  of  the  terms  [xloowp  ■  (t)  *0  (■y)  Is  depend  upon  solar 

radiation,  air  temperature^  and  vapor  pressure  inoresients. 


. . 


(4.)  "  total  hemiapherieal  radiation  per  unit  area  at  the  ground 

X/  ♦.Vj 

^  ^  "  ground  radiosity  per  unit  area 

tq  I  -  heat  transfer  rate  per  unit  area  into  or  out  of  the  ground 
A  fig  at  the  interface  ■ 

' '  I  ^  1  "  transftf*  rate  per  unit  area  due  to  evaporation  or  coo 

I  A  /  ^  denaatlon  of  water  vapor  at  the  interface 

la  the  abaence  of  direct  measurenenta  of  either  the  eonveotive  or  latent 
2 

heat  tranafer  ratea  the  following  method  of  aeparating  the  aum  of  theae  two 
terma  ia  auggeatad.  Equate  the  defining  convective  and  latent  heat  tranafer 
equationa  within  the  turbulent  lagrer  adjacent  to  the  laminar  aublagrer  to  the 
aum  under  conaideratioa. 


"  ground  radioaity  per  unit  area 


heat  traaater  rate  per  unit  area  into  or  out  of  the  ground 
at  the  interface- 

heat  tranaftf*  rate  per  unit  area  due  to  evaporation  or  con- 
denaatlon  of  water  vapor  at  the  interface 


+  A. 


•  -ycj.*. 


ax  + 


vAiere* 


>  vertical  eddjr  diffuaivitjr  for  aaaa  tranafer 
<1^  a  vertical  eddy  diffuaivity  for  heat  tranafer 
C  a  water  vapor  concentration 

1  a  latent  heat  of  vaqpMleatioQ 

Several  inveatigatora  in  the  literature  have  ahown  that  the  heat,  maae,  and 
aommtum  eddy  diffuaivltiea  are  deadly  related  (Refezwneee  9  nad  10).  It  i« 
thua  propoaed  that  the  heat  and  maaa  tranafer  eddy  diffusivitiea  are  related  to 
each  other  by  a  eonatant,  j.  That  la. 


^apafimeterr  of  the  type  uaed  by  Ibomthwaita  (Refemee  8)  would  appear  to 
be  r0tf  uaefhl  ia  meeavaiag  mean  araporatloa  ratea  if  not  inatantaneoua  onaa. 
Kaaa  flLow  aMtera  aii^ar  la  prtaol^e  to  heat  flow  aetera  are  eurrently  being 
oonaidered  by  the  Ohlveraity  of  Oalifomia,  loa  JUg^aa* 


> 
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'  Upon  the  substitution  of  equation  (8)  into  equation  (7)»  the  heat  eddy  diffusivity 
can  be  expressed  as 


References  9  and  10  indicate  that  for  duct  flew  systems,  the  constant,  J,  is 
equal  to  unity.  Thus  an  approximate  method  of  separating  the  convective  and 
evaporative  heat  floir  terms  consists  of  postulating  that  the  constant  j  ■  1, 
solving  for  by  equation (9),  and  hence  evaluating  the  convective  and  evapora¬ 
tive  heat  flow  terms.  The  ratio,  j,  will  become  more  firmly  established  with 
continuing  micrometeorological  research.  In  the  event  that  one  of  the  two  heat 
flow  terms  is  small  compared  to  the  other,  no  separation  problem  exists. 

After  having  determined  the  convective  heat  transfer  rate,  it  is  possible 
to  determine  the  eonvect.ive  neat  transfei-  uuitJuetaneea  and  thermal  eddy  diffu- 
sivities  from  the  defining  basic  eddy  diffusion' relations  which  have  been 
presented  previously.  Itass  transfer  conductances  and  dlffusivities  may  be  ob¬ 
tained  In  a  similar  manner. 

Some  limited  convective  heat  transfer  measurements  and  calculated  thermal 
conductances  are  presented  in  a  following  section. 

B.  Vertical  Momentum  Transfer 

Boundary  drag  foirce  measurements  in  pipe  flow  systems  can  be  made  with 
relative  ease  in  comparison  to  boundary  drag  force  measurements  in  the  atmos¬ 
pheric  system.  No  direct  measurements  of  the  drag  force  exerted  hf  the  wind 
on  the  earth's  surx'ace  appear  to  be  reported  in  the  literature  with  the  excep¬ 
tion  of  those  of  P.A.  Sheppard  (Reference  11). 

The  preblen  of  measuring  the  air  shear  stress  at  ^he  earth's  surface 
essentially  consists  of  determining  Uie  dtsg  force,  on  a  small  area  of  that 
surface  1)  whose  surface  characteristics  are  r^>resentatlve  of  the  surrounding 
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surface  and  2)  which  Is  so  located  that  the  velocity  profile  above  it  is  typical 


of  the  surrounding  velocity  profile  (no  new  boundary  layers  initiated).  The 
atmospheric  shear  meter  that  was  developed  by  the  University  of  California  was 
designed  with  the  above  two  requisites  in  H-'ni  s-id  has  yielded  encouraging 
results.  This  shear  meter  essentially  .:onslsted  of  a  large  tank  filled  with 
water  into  which  was  placed  a  shallow  float  whose  surface  was  representative 
of  the  surroundings.  Drag  forces  were  measured  by  noting  the  defle<'tion  of  a 
sensitive  coil  spring,  one  end  of  which  was  attached  to  the  float  ard  the  other 
end  to  the  rim  of  the  tank.  Abrupt  flow  discontinuities  were  avoided  by  filling 
the  latter  to  within  less  than  l/l6  of  an  inch  of  the  brim  with  water  and  adjust¬ 
ing  the  .weight  of  the  float  so  that  it  extended  less  than  lA^  inch  above  the 
water  surface.  The  tank  rim,  float  surface,  water  surface,  and  surrounding 
earth  surface  were  thereby  essentially  in  one  plane.  Smoke  flow  studies  indi¬ 
cated  that  smooth  air  flow  conditions  existed  over  the-shoar  motorr - 

The  shear  meter  construction  details  are  shown  in  Figure  1.  The  tank  was 
U8  Inches  in  diameter  and  3  inches  deep;  the  respective  dimensions  of  the  float 
were  36  and  2  l/U  inches.  It  was  necessary  to  cover  the  upper  float  surface 
with  a  thin  layer  of  earth  in  order  to  create  a  sxurface  which  was  similar  to  the 
surrovindings .  fhe  total  weight  of  this  float  was  about  seventy-five  pouads. 
Because  it  was  desired  to  measure  the  mean  shear  stress  rather  than  the  instan¬ 
taneous  values,  a  coil  spring  with  a  very  low  spring  constant  was  utilized  in 
order  that  the  shear  fluctuations  superposed  upon  the  mean  shear  would  not  causa 
the  float  to  oscillate.  The  coil  spring  was  made  by  winding  piano  wire  (0.007 
inch  diameter)  on  a  steel  rod.  One  end  of  the  spring  was  attached  to  the  rim 
of  the  ta{^  by  a  clip  idilch  could  be  rotated  idien  significant  changes  in  wj.nd 
direction  occurred.  The  other  end  of  the  spring  was  fixed  to  a  wire  located 
under  the  float  and  attached  to  the  float  center.  This  arrangement  allowed 
angular  freedom  of  the  float  without  producing  spring  ext«islon  due  to  angular 


•  t 
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float  rotation. 


I 

1 

The  drag  force  for  this  simple  spring  system  is  given  l?jr 


where, 

-  spring  constant 
8^  -  deflection  of  the  spring 

\  -  shear  meter  surface  area  exposed  to  air  flow 

The  spring  constant  for  the  particular  spring  used  was  0.001  Ibs/inch  and  was 

a 

accurate  within  one  percent. 

Although  the  shear  meter  that  has  been  di..cribed  was  devdoped  ind^ai- 
dently  of  Sheppard's  meter,  the  two  instruments  are  similar.  However,  Sheppard's 
meter  was  based  on  a  torsion  principle  rather  than  a  tension  principle, and  also 
his  meter  was  designed  to  measure  drag  forces  over  smooth  surfaces.  Field  exper¬ 
ience  with  the  authors'  present  shear  meter  has  Indicated  that  the  following 
la^rovements  are  desirable!  1)  a  response  system  which  is  independent  of  wind 
^Ire^ion  and  2)  a  shear  stress  recording  mechanism.  Development  of  future 
instruments  having  these  in^>rovements  has  been  initiated. 

After  having  determined  the  shear  stress  with  the  aid  of  a  shear  meter, 
it  is  possible  to  calculate  the  drag  coefficient  and  momentum  eddy  diffusivity 
from  the  defining  eddy  diffusion  relations  idiich  have  been  presented  previously. 

Some  lladted  shear  stress  measurements  and  ealcxilated  drag  coefficients 
and  eddy  diffusivitles  are  presented  in  a  following  section. 


V. 


m 


ANALYTICAL  AND  MUKERICAL  EDDY  DIFFUSIOW  ANALYSES 
A.  Analytical  Analyses 

A  nuaber  of  analytical  eddy  diffusion  analyses  can  be  found  in  the  litera¬ 
ture.  Prandtl  (Reference  12)  derived  the  steady  state  logarithmic  velocity  pro¬ 
file  for  the  surface  layer.^  Rossby  and  Montgomery  (Reference  13)  and  Sverdrup 
(Reference  14)  developed  steady  state  tdnd  velocity  profiles  for  the  surface  layer 
for  adiabatic  and  stable  atmospheres.  Letbau  (Reference  15)  has  recently  derived 
steady  state  velocity  and  temperature  profile  expressions  for  the  surface  layer 
for  adiabatic  and  non-adlabatic  atmospheres.  0.  0.  Sutton  (Reference  16)  has  devel¬ 
oped  point  and  line-source  eddy  diffusion  solutions  for  stable  atmospheres.  Brunt 
(Reference  17)  has  presented  a  periodic  convective  heat  flow  solution  for  a  system 
in  which  the  boundary  temperature  varies  sinusoidally  with  time,  thus  approximating 
the  diurnal  soil-air  Interface  temperature  variation;  the  eddy  diffusivlty  was 
postulated  to  be  Independent  of  time  and  height.'  Haurwitz  (Reference  IS)  has  devel¬ 
oped  a  periodic  convective  heat  flow  solution  for  a  system  in  which  the  boundary 
temperature  varies  sinusoidally  with  time  and  the  eddy  diffusivlty  varies  linearly 
with  height  but  is  Independent  of  time; 

In  the  following  paragraphs,  two  new  periodic  convection  analyses  in  the  ab¬ 
sence  of  atmospheric  thermals  are  presented.  One  analysis  pertains  to  an  eddy 
diffusion  system  in  which  the  boundary  temperature  varies  sinusoidally  with  time 
and  the  eddy  diffusivlty  varies  sinusoidally  with  time  (to  approximate  .unstable 
diffusion  during  the  day  and  stable  diffusion  at  night.)  Another  analysis  pertains 
to  an  eddy  diffusion  system  la  which  the  boundary  temperature  varies  sinusoidally 
with  time  and  the  eddy  diffusivlty  varies  sinusoidally  with  time  a.)d  linearly  with 
height. 

1)  Eddy  Diffusivlty  Varies  Sinusoidally  With  Time  But  Not  With  Height  < 

The  heat  transfer  differential  and  boundary  erpiations  for  the  system  under 

3  ' 

The  surface  layer  has  been  defined  as  the  air  layer  adjacent  to  the  earth's 
surface  where  the  vertical  convective  heat  flow  and  vertical  fluid  shear  stress 
are  uniform  with  height. 


eon sidsra lion  are^ 
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where. 


3  t  ,  3 

(e 

(11) 

3d  3  s 

- 

t(0,  6)  » 

tg  MS  Oi  6 

(12) 

t(«.  «) 

m  0  • 

(13) 

S  ‘  ‘•i 

+  Cg  eos  ft»  6 

(14) 

c  >>  potential  temperature 

■  amplitude  of  the  sinusoidal  boundary  temperature  ware 

<•>  ••  frequency  of  the  sinusoidal  boundary  temperature  and 

dlffuslvlty  waves 

^  ‘  period  of  the  sinusoidal  boundary  temperature  and 
dlffuslvlty  waves 

cj.  «  constants  In  equation (14) 

If  the  differentiation  Indicated  In  equation  (11)  Is  performed,  the  following 
equation  results: 


1  3t  .  /  l 

,  (15) 

3  6  +  CgeototSjZO 

3*» 

Equation (15) can  be  simplified  by  making 

a  change  of 

variable. 

Let 

■ 

-^$ino>9 

<a 

(16) 

Then 

-  c,  ♦ 
3d  ^ 

CM  0)6 

(17) 

Now  the  left  hand  term  of  equation (15) can  be  simplified  as  follows: 
I  I  \  3 1  .  I  1  )  3A  Lt  - 

+  Cj  oof  a 6/  ^6  "i  K _ 3 \ _ 


(W) 


^ese  equations  are  expressed  in  texmis  of  the  potential  temperature  rather  than 
the  ordinary  temperature.  Potential  teaqMroture  ean  be  expressed  ar  ^  7  *  rs. 
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The‘  transforaed  dl^erentlal  and  boundary  equations  are 


t  (0,  X.)  «  t-  eoa  ct>  /(X) 


-  0 


(30)  ' 


iihere«  y(X)  ■  6 

It  can  be  sho'm  that  the  function  c  (0.  X)  is  an  even  periodic  function  with  a  period 
00  Tvhere  0^  is  the  period  of  the  sinusoidal  boundary  temperature  and  diffu- 
sivity  variations.  The  function  t  (0.  X)  nay  thus  be  expressed  as  the  Fourier  cosine 


series 


where, 


t  (0,  X) 


a 

Ar' 


(X')  dX' 


2  t  (X.')  dX' 


Equation  (19)  can  be  solved  by  the  separation  of  variables techniq:ae. 

( 

Let  t  ■  Z(x)A(X)  (23) 

'where  Z(x)  andA(X)  are  functions  of  x  and  X,  respectively.  The  steady  state 
periodic  solution  of  equation  (2?)  which  satisfies  ecpatlons  (20)  «nd  (21)  is 

Z(x)  ■  eon«t«nt|  ,  A(X)  ■  opnatwvS^ 

or  s  -  (20 

^ . r . .  . - 


i.*-  — r 


Id 


The  transformed  differential  and  boundary  equations  are 


3 1  ,  3»  e 

^  77  . 

(19) 

t  (0,  X)  -  oos  ci  t(S) 

(30) 

,*i*,  t  (a.  X)  -  0 

(21) 

where,  /^)  ■  0 


It  can  be  shown  that  the  function  t  (0,  X.)  is  an  even  periodic  function  with  a  period 
0^  where  0^  is  the  period  of  the  sinusoidal  boundaiy  temperature  and  diffu- 
slvlty  variations.  The  function  t  (O.  X)  may  thus  be  expressed  as  the  Fourier  cosine 
series 


(0,  X)  .  ^ 


(22) 


where. 


*0  ‘  ^  I  ^  t(X')dX' 


o/a 


oo« 


2  nTT  V 


aX' 


Equation  (19)  can  be  solved  by  the  separation  of  variables  technique. 


Let  t  ■  Z(z)A(X)  (23) 

'where  Z(t)  andA(X)  are  functions  of  a  and  X,  respectively.  The  steady  state 
periodic  solution  of  equation  (23)  which  satisfies  equations  (20)  and  (2l)^ls 
Z(a)  •  coMtonSj  ,  A(X)  ■  ^IvX 


O? 


_ _  _ _ _ _ .A..  *•. 

«h«re«  0  and  v  ar«  constants  which  arc  to  be  dstsnnined  and  I  ••  ,  Both 

tha  ra^  and  laaginarj  parts  of  the  convex  solution  given  in  aquation  (24)  ara 
solutions  of  tha  original  diffarantial  aquation*  On  physical  grounds,  tha  real 
part  of  aquation  (24)  unist  ba  chosan.  Also,  tiia  solution  must  ba  axprassad  in 
tanas  of  a  serias  baoausa  this  is  tha  form  of  tha  boundary  aquation  (22). 


Ihua,  t(i,  X)  ■  Ci,  •  ^  ^  ”  i?  *) 


whan  n  >  0,  1,  2, 


I  Tha  eonstants  and  ara  datarminad  from  tha  botnadsry  .aquation .  (22) . 


At  a  ■  0. 


K)  -  ^  J  ••  ••*  •  S  *  ••*  **% 

'  an  *1  %  a^i - 


Thus, 


C.  •  3 

^  2 


•  %  for  n  >  I,  2, 


Tha  tanparatora  solution  for  tha  hast  transfar  aystam  undar  eonsidaraticn  is 

T(s,  0)  •  T  a  ♦  a  sfo  a| 

whara,  a  ■  0,  1,  2,  •*• 


,\  . 


~L  I  2  t(\')  CM  2n^' 


The  periodic  vertical  heat  flow  solution  for  the  heat  transfer  system  under  considera¬ 
tion  can  be  obtained  by  substituting  the  temperature  solution  givmi  in  equation  (27) 


into  equation  (1) 


.  -yc.elll  *  r) 

^  A  jeonp  r  *f\  a  X  / 


Mass  concentration  and  mass  flow  rate  solutions  Can  be  derived  Just  as  the  tem¬ 
perature  and  heat  flow  solutions  were  drived  above  if  the  mass  transfer  system 
possesses  the  above  stipulated  boundary  and  eddy  diffuslvity  equations. 

2)  Eddy  Dlffhsivity  Varies  Sinusoidally  With  Time  end  Linearly  With  Height 


The  heat  transfer  differential  and  boundary  equations  under  consideration  are 


3  3t\ 


t(0. 6)  -  t.  eo*  <o3 


^  ■=  0 


where. 


■  (bj  +  *)  (l  +  hg  eo*ci>0)  (31 

bj,  b^,  and  b^  are  constants  in  equation  (31) 

If  the  differentiation  indicated  in  eqxiation  (28)  is  performed,  there  results 
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the  following  equation 


3<- 


3t  .  t  «  9Ls 

Td  -ST  n 


•  (1  +  bj  eos  <0  6)  ^  +  (b^  +  b^  x)  (1  ♦  hg  ee*  CJ$) 


or 


■  >'.1;  ' 

Equation  (32)  can  be  simplified  by  making  a  change  of  variable. 


Let 

That 


&  +  —^sindO 


1  60*  <t>  B 


The  left  hand  term  of  equation  (32)  cun  be  simplified  as  followst 


(32) 

(33) 

(34) 


(35) 


The  transformed  differential  and  boundary  equations  are 


U  ■ 

(36) 

t(0.  ,3)  -  C05  w  /  (|S) 

(37) 

o 

■ 

II 

M 

(38) 

where*  f{^  •  B 

The  function  t(0,  fl)  Is  an  even  periodic  function  with  a  period 

m«y  be  expressed  as  a  Fourier  cosine  series* 

^  and  thus 

t(0.  ^  ^  eo* 

(39) 
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Sir 


•where. 


go  -  2  t(/3')  ay? 

0-& 


I 


f 

r-» 
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A  further  change  of  variable  irlll  simplify  the  differential  equation  (36). 


Let 


Then  equation  (36)  reduces  to 


**1  +  *»,  « 


3 1  ,  ^  ^  t  it  Ls 

S  /J  ^  3  t*  4  r  3  e 


(40) 


(41) 


TTguatlnri  (1|1  )  t»*n  tiw  Hnivnrf  hy  tha  naparnt.inn  of  variable« technique.  _ 

Let  ,  t  -  R(c)  B(fl)  (42) 

where  R(r)  and  B(/3)  are  functions  of  t  and  respectivdy.  The  steadjr 
state  periodic  solution  of  equation  (Ul)  idtich  satisfies  equations  (37)  and  (38) 
is 


R(r)  •>  constant^  Ji  cj. 


m 


conaUmtg  P  ^ 


or 


t  - 


(43) 


where  ^  and  P  are  constants  to  be  determined  and  where  the  function 
1^-U  ji  can  be  expressed  as 


-a[.«(^r)  *  or 

(Reference  19).  The  tera  is  the  modulus  «id  the  term 


airplitude  of  the  eoa:plex  function 


p.). 


> 


That  is. 
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and 


tan  < 


Thus  equation  (U3)  ma^  be  expressed  as 


k«r 


ktl 


(f'l 

w 


(44) 

(45) 


in  order  that  the  boundary  condition  car.  be  satisfied,  that  is 


t(»  -  0) 


(46) 


must 


be  of  the  form  «"^P^  (the  real  part  being  eoa  p  ). 

t  b, .) .  -p<j 

where,  is  a  constant 

The  real  part  of  solution  (U?)  is  chosen  on  physical  grounds.  Also,  the  solution 
must  be  expressed  in  tezms  of  a  series  as  is  the  boundary  equation. 

Thus, 
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B,  .  l.'urisrical  Analysea 

Several  numerical  methods  of  evaluating  eddy  diffusivity  profiles  from  tran¬ 
sient  air  tenperaturo  and  humidity  profiles  in  the  atinc sphere  Lave  been  given  in 
the  literatiire  (RejTerences  20  and  21,  for  exaripleN  These  methods  usually  involve 
replacing  the  differential  equation  describing  transient  diffusion  by  a  finite 
difference  equation.  HThen  experimental  temperature  or  humidity  data  are  substituted 
into  the  finite  difference  equation,  a  solution  can  be  effected.  Some  of  these 
methods  involve  procedures  which  do  not  appear  to  be  sufficiently  general.  For 
exanple,  one  method  requires  that  at  one  time  in  the  analysis,  eddy  diffusivlties 
in  adjacent  air  layers  be  the  same.  Another  method  requires  an  error  analysis  to 
be  applied  to  the  resulting  diffusivity  data  because  of  an  over  simplified  postulate 
that  was  used}  this  error  analysis  then  invalidates  much  of  the  data. 

- Thrpa  mimerical  methods  of  datemining  eddy  diffusivlties  as  a  function  of 

height  and  time  from  transient  temperature  and  humidity  profiles  are  presented 
below.  These  methods,  which  do  not  appear  to  be  presented  in  the  literature, 
possess  several  desirable  advantages. 

l!ethod  Ko.  1 

A  heat  rate  balance  on  a  vol\ime  of  air  -  Zj  units  thick  with  a  unit  area 
base  is 


*  '’1 


ds 


(51) 


It  is  new  postulated  that  the  eddy  diffusivity  varies  linearly  with  height,  that  le 


-0  ♦  (52) 

where  and  are  constants  for  a  given  time.  If  it  is  desired  to  determine 

the  eddy  diffusivity  profile  (a  linear  one  wliich  approx^ates  the  actual  one)  in  a 
given  air  layer*  the  following  procedure  is  suggested.  Divide  the  laser  in  question 
into  two  layers  and  write  equation  ($1)  for  each  of  them.  Upon  substituting  the 
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linear  diffuslTlty  rd.ation  given  by  equation  ($2)  into  the  two  heat  balance  equa- 
tlona, there  reault  the  following  two  equations t 


».('•  *  tXH  *  ''1  ■  ’i  ‘‘KH  ’  '"I 


If  density#  heat  capacity*  adiabatic  l^se  rate*  and  teaperature  gradient  data  for 
a  given  set  of  experinental  temperature  measurements  at  a  given  time  are  substituted 
into  equations  ($3)  and  (^U)*  and  if  the  int^rations  indicated  are  performed*  there 
result  two  algebraic  equations  in  two  uidcnowns  (c^  and  «j).  The  solutions  of 
theae  equations  are  eaeUy-obtained^ 


Method  No.  2 

k  somewhat  more  general  eddy  diffuslvity  distribution  can  be  e3q>ressed  by  the 
following  power  ejqpressiont  • 


«  •  «-  * 


‘e  • 

where  e^  and  are  constants  for  a  given  time.  Upon  substituting  equation 

{$$)  into  the  two  heat  balance  equations  and  inserting  the  experiments,  property 
and  temperature  data*  as  was  indicated  in  Method  No.  1  dbove*  two  solvable  Sge- 
bralo  equations  in  two  unknowns  result. 

Msthod  Nd.  3 

in  eddy  diffiielvlty  distribution  idtioh  is  very  general  is  the  following  series 

s  *  ......  - - 

eontalning  hi*  tenHi 


<  >  V  *  *1*  *  *  ••• 


(*) 


V  *j*  •••  *»-i 


are  eonstants  for  a  given  time,  dividing  the 
air  layer  under  oonftderation  into  *tt'  email  er  ligrers  and  substituting  eqnatien  ..r 
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C?6)  Into  the  'n»  heat  balance  equations. 
In  unknowns,  if  .n.  i,  a  large  number 
the  aid  of  computing  machines. 


there  would  result  'n*  algebraic  equations 
,  these  equations  could  be  solved  with 


Methods  1,  2,  and  3  are  thus  sug/jested  as  ♦  u  . 

Buggesued  as  rapid  techniques  of  determining 

transient  eddy  dif fusivity  profiles  for  ^ 

y  promts  for  heat  and  mass  transfer  from  transient 

temperature  and  moisture  data  ^ 

data.  These  methods  are  founded  on  general  postulates. 


AHALYSI3  OF  SEVERAL  SETS  OP  IHCRCgfflTEOROLOQICAL  DATA 


A.  Experimental  Determination  of  ConveetlTe  Heat  Flows*  Convective  Conduotanoes* 
Shear  Stresses*  Drag  Coefficients*  and  Eddy  Diffusivities* 

1.  Heat  Transfer 

On  a  clear  summer  day  between  the  hours  of  lOtOO  A.1I.  and  3tOO  P.H .  on  August  12 
I9U9*  a  preliminary  air-earth  Interface  heat  balance  study  was  conducted  nepr  Tan 
Nuys*  California  (Reference  22) .  This  study  involved  the  measurement  of  1)  air  tem¬ 
perature  profiles  in  the  lower  four  feet*  2)  earth  temperature  profiles  in  the  upper 
ten  inches*  3)  solar  irradiation*  U)  ground  heat  flows*  and  $)  wind  velocities  At 'A 
five  foot  height.  The  eiqperiiKental  site  was  located  on  a  large  flat  field  covered 
with  scattered*  short  brush  less  than  a  foot  hi|^.  Thermocouples*  a  directional 
radiometer  (Reference  23) «  a  heat  meter*  and  a  cup  anemometer  were  used  to  meAsure 


ieraiure* 


Locity*  respectively. 


The  interface  heat  balance  relation*  equation  (6)  *  may  be  expressed  as 


'  X 


(x)<i  -  “'(jj-i-  ’  ■  ‘V"''*  *  (ij>, 


The  total  hemispherioal  radiation  and  ground  rediosity  have  been  expressed,  in  terms’ 
of  solar*  gaseous*  and  ground  radiation^.  The  ayabols  a  nd  are  the 

grey-body  eax^h  surface  emissivity*  Rtefan-Boltssann  constant*  and  the  interface 
tanperatore*  respectively.  Because  the  ground  was  diy  and  had  not  been  exposed  to 
rain  for  more  than  a  two  month  period  the  evaporative  heat  losa  term|^|^  was 
postulated  to  be  small  compared  to  the  other  texiw  In  the  equation.  The  solar  rad¬ 
iation  and  ground  host  flow  terns  wore  measured  by  the  Instruments  described  abort'* 
Data  In  the  literature  and  some  recant  nsasiO'SBents  by  the  Tlmnaal  Badlatlcn  Frojast 
at  tile  Q^reralty  of  CAUfomla  Indicated  that  ^  was  very  close  to  0.6.  The 

tern  I  ~  ?  was  estlncted  ftrcnnoetumsl  (gaaeoM) 


^ad  a  heniqphericAl  radleneter  rather  than  a  dlreetionsl  rarticawter  bean  araildULt 
at  tiM  expsaimentai  site*  this  re-arrsngenent  of  equation  (6)  would  not  hare  been 
nioessaiy*-:' .  -  . _ ; _ ;  _ 
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radiation  calculations  (Reference  2lt)  and  the  interface  temperature  measurements.  This 
estimate  was  in  agreement  with  measurements  by  Brooks  and  Kelly  (Reference  20),  It 
was  thus  possible  to  calculate  the  convective  heat  loss  using  equation  (07) » 

The  results  of  the  Van  Nuys  heat  transfer  study  are  presented  in  Figiires  2, 

3.  U,  and  0.  The  unit  thermal  conductance,  calculated  by  equation  (l»)  and  shown 
in  Figure  2,  is  observed  to  be  approximately  proportional  to  the  wind  velocity 
(Figure  3).  It  is  interesting  to  recall  that  the  thermal  conductance  for  turbulent 
flow  in  ducts  and  over  flat  plates  varies  as  the  eight  tenths  power  of  the  fluid 
velocity.  The  magnitude  of  the  thermal  conductances  were  in  agreement  with  certain 
measurements  by  F,  A.  Brooks  (Reference  20),  which  v/ere  made  under  similar  circum¬ 
stances.  Figure  U  reveals  the  variation  of  air  temperature  difference  between  the 
surface  and  foiir  feet  with  time.  Figure  0  shows  the  vertical  convective  heat  flow 
which  appears  to  have  a  maximum  near  one  o'clock;  this  maximum  almost  coincides 

with  the  solar  irradiation  maxinmm  (at  12t30  P.M.)  and  the  surfaca-air  temperature - 

increment  maximum  (at  12:30  P.M.). 

2.  Momentum  Transfer 

Preliminary  momenttun  transfer  studies  in  the  lower  twenty  feet  of  the  atmos¬ 
phere  were  conducted  at  intervals  during  January  and  February,  19^0,  at  Riverside, 
California,  These  studies  involved  the  measurement  of  wind  velocity  profiles  and 
atmospheric  shear  stresses.  Sensitive  cup  anemometers  and  the  shear  meter,  which 
has  previously  been  described,  were  used  to  measure  wind  velocity  and  shear  stress, 
respectively.  The  experimental  site  was  located  on  a  flat  rectangular  field  which 
had  a  disked  surface,  settled  and  flattened  by  rain.  The  few  rough  portions  of 
ground  that  did  exist  within  ten  to  twenty  feet  of  tlie  shear  meter  were  smoothed 
with  a  rake.  Ro  vegetation  covered  this  field,  although  some  scattered  patches  of 
grass  wore  noted  on  the  last  day  of  operation  (February  20th) .  Most  of  the  experi¬ 
mental  measurements  were  made  in  afternoon,  between  2:00  and  lit 30;  atmospheric 
conditions  of  slight  instability  existed  at  these  times.  A  few  measurements  were 
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f  6 

I'  also  made  under  density  flow  conditions”  between  6it(U  A.M.  and  6:^  A.K.  Figure  6 
^  shows  a  graph  of  shear  stress  as  a  function  of  wind  velocity  at  the  eighty  inch 

I  level.  The  shear  stress  and  velocity  values  were  averaged  over  one  to  five  minute 

-  r 

f  time  intervals.  The  afternoon  shear  stress  data  were  noted  to  vary  as  the  square 
of  the  velocity  and  were  represented  by  equation  (^);  the  drag  coefficient  was  found 

: 

r  to  be  0.00137 •  This  relation^  which  is  noted  as  curve  ’'A"  in  Figure  6,  represents 
the  data  within  an  accuracy  of  ♦  309S.  Unfortunately  the  morning  data  consists  of 
only  two  points.  These  points,  however,  are  believed  to  be  more  accurate  than  some 
of  the  afternoon  points  because  the  stable  air  flow  conditions  at  that  time  yielded 
very  steady  shear  stress  and  wind  velocity  measurements.  A  square  velocity  relation 
faired  through  these  two  points  and  the  origin  is  shown  by  curve  "B"  in  Figure  6. 

It  is  intended  that  this  c\rrve  suggest  the  effect  of  stability  on  the  shear  stress  - 
wind  velocity  relation)  many  sets  of  data  must  be  obtained  before  the  stability 
effect  can  be  thoroughly  studied.  The  exponent  on  the  common  power  wind  velocity 
relation  is  about  0.16  for  curve  "A"  and  about  0.33  for  curve  "B".  A  calculation 
has  revealed  that  If  it  were  possible  for  conpletely  laminar  flow  to  exist  tuider 
the  existing  velocity  conditions,  a  shear-velocity  curve  falling  far  below  curve 
"B”  would  result.  If  further  research  substantiates  the  above  trend,  it  appears 
that  for  a  given  wind  velocity,  the  atmospheric  shear  stress  at  the  ground  would 
Increase  with  a  decrease  in  the  power  law  exponent  (increased  turbulence).  This 
behavior  would  seem  reasonable  in  the  light  of  our  present  knowledge  of  fluid 
turbulaice. 

A  comment  on  the  height  at  which  the  wind  velocity  is  to  be  measured  and 
the  accompanying  implications  seems  to  be  in  order.  Consider  the  two  velocity 

■ 

profiles  shown  in  Figure  7}  one  profile  is  stable  and  one  is  unstable.  Also  con- 
-sider  that  the  stable  profile  is  associated  with  a  smaller  shear  stress  than  is 
the  unstable  profile.  _  _ 

^loctumal  cooling  of  the  earth's  surface  creates  tenq>erature  inversions.  An 
atmospheric  system  with  such  inversion  conditions  and  sloping  terrain  then 
possesses  ai  density  potential  which  causes  the  cool  air  in  the  air  layers  con- 
tiguooe  ta^e  gpround  to  now  fton  higher  to  lower  elevations. 


Sktar  Stress]  tlft^ 
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Figure  6*  Sheer  Stress  es  «  Funotion  of  - 

Wind  Velocity  at  SO  inches.  ^  | 

- . _ _ i 

^  f 
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,  Figure  7*  Stable  and  Unstable  Wind  Velocity  Profiles 

Consider  the  shear-velocity  relations  that  result  for  each  of  the  tiro  proposed 
helots,  Zj  and  zg  for  the  two  velocity  profiles  shown.  If  Zg  Is 
chosen  as  the  height  at  which  the  velocity  is  to  be  measured,  it  is  seen  that  the 
stable  shear  stress  is  less  than  the  unstable  shear  stress  for  the  sane  wind 
velocity.  However,  if  some  lower  reference  height  such  as  Zj  is  chosen,  the 
stable  velocity  is  so  much  less  than  the  unstable  one  that  the  stable  and  unstable 
shear-velocity  relations  now  may  be  much  closer  to  each  other  than  they  were  when 
the  reference  zg  was  used.  Further,  if  one  chooses' the- reference  height  so' that 
it  lies  within  the  laminar  sub-layer,  a  single,  linear  shear-velocity  curve  would 
exist  for  stable,  neutral,  and  unstable  flow.  That  is,  within  the  laminar  sub-layer, 
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A  ooiiqparlaon  of  the  Riverside  shear  data  with  some  of  the  shear  datA, 
reported  in  the  literature  is  made  in  Figure  8.  All  velocities  have  been  re¬ 
ferred  to  a  reference  height  of  thirty  feet.  In  some  eases  it  was  necessary  to 
determine  the  velocity  for  the  new  refermce  height  when  complete  velocity  data 
were  not  reported}  the  seventh  power  velocity  relation  was  utilised  to  make  the 
determination. 

I 

Sutcliffe  and  Taylor  have  utilised  the  terhnioue  of  evaluating  boundary  « 

shears  from  the  hydrodynamic  equations  that  describe  the  vertical  wind  spiral  up  ; 

to  the  gradient  wind  level}  experimental  vertical  wind  velocity  measuresisnts  are 
neoessary  to  make  this  determination.  This  method  is  ]>3redieated  on  the  postti-- 
lates  that  steady  state  flow  conditions  are  established  and  that  unidirectional 
flow  exists  at  the  earth's  surface.  Shear  stresses  calculated  by  this  method  rep¬ 
resent  mean  values  for  the  earth's  surface  and  include  the  drag  effects  of  such _ 

regiona  as  bodies  of  watery  plowed  fields,  brush  land,  orchards,  and  foreate. 

Suteliffo'e  daita,  curve  "A*,  and  Taylor's  data,  curve  "B*,  which  represent  shear 
data  over  the  Salisbury  Plain,  fall  within  forty  percent  of  one  another* 

Sutcliffe  also  has  obtained  Aear  data  over  the*  South  Atlantic  with  the  aid 
of  kites  (see  ourve  "B"}}  it  seems  reasonable  that  this  ourve  should  fall  far 
below  curves  "A"  and  "B"  because  the  ocean  is  a  relatively  smooth  surface  oonparod 
to  land* 


Curve  "C  represents  a  shear-velocity  relation  which  Rossby  obtained  by  sub¬ 
stituting  ejqwrinsntal  wind  velocity  data  over  Uie  ooean  into  von  Karmm's  goner- 
lUsad  vslooity  dietrlbutiSR  for  flow  In  a  wooth  duet  syirtem. 

The  direct  shear  measurements  of  Sheppard  with  his  Shear  meter  are  noted 
by  (^irvs  These 


surements  were  mads  over  a 


th  oonerete  surfaee* 


It  is  of  interest  be^onpare  Sheppard**  resslts  with  these  of  the  aathors*; 


as  the  flow  systmas  were  somsuhat  similar  and  bsoauss  direct  methods  of 
shear  stims  wmi«  used.  It  is  suggested  that  Sheppard's  data  lie  above  the 


authere*  data  Wsause  of  the  diff ersnee*  in  atmsspherie  etahlUty  in  tho  two  eases," 
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1 

The  velocity  pbwer  lavr  exponent  for  Sheppard's  data  was  about  0.07  lyhereas  the 
exponent  for  the  author's  data  was  about  0.16 » 

Momentum  eddy  diffuslvitles  have  been  determined  from  the  Riverside  data  with 

i 

the  aid  of  equation  (2).  The  three  typical  profiles  shown  in  Figure  9  indicate  that 
the  eddy  diffusivity  la  approximately  proportional  to  the  wind  velocity.  Note  that 
one  can  show  from  Prandtl's  analytical,  steady  state,  logarithmic,  velocity  solution 
that  the  eddy  diffusivity  is  directly  proportional  to  the  velocity.  Note  also  that 
these  dlffuFivity  profiles  have  configurations  similar  to  those,  existing  in  the  vlo~ 
inlty  of  the  wall  in  duet  flow  systems;  the  diffusivity  varies  almost  linearly  with 
distance  from  a  very  small  value  (molecular  viscosity)  at  the  boundary. 

B.  Numerical  Eddy  Diffusion  Analyses  of  Transient  Air  Temperature  Data 

Thermal  eddy  diffusivity  profiles  have  been  evaluated  for  two  sets  of  transient 
air  temperature  measurements  by  means  of  numerical  method  No.  1  which  has  been  des> 
crlbed  above.  One  set  of  data  waa-0btalaed-4uring-a-TitffaBlen~lHwstigatlon  in  the 
desert  at  Datelan,  Arizona,  (Reference  29)  and  another  set  was  obtained  on  gentlyr 
rolling  cotton  land  at  Manor,  Texas,  (Reference  30). 

Figure  10  shows  some  typical  thermal  eddy  diffusivity  profiles  as  a  function  of 
time  for  Datelan,  Arizona.  Calculations  were  not  made  for  the  late  morning,  noon', 
and  early  afternoon  periods  because  atmospheric  thermals  existed  at  these  times^. 

The  diffusivity  profiles  in  Figure  10  indicate  the  influence  of  the  diurnal  turbul¬ 
ence  variation  on  convective  heat  transfer  (high  diffuslvitles  dwing  the  day  and 
low  ones  at  night).  The  diffusivity  profile  at  0450  has  been  further  sub^jivided 
(note  broken  portion  of  curve)  in  order  to  show  the  influence  of  noctximal  stability 
on  the  diffusivity.  The  presence  of  this  stable  region  was  first  suggested  by  the 
extreme  temperature  inversion  conditions  above  the  one  hundred  foot  level.  Note  that 
all  diffusivity  pz*ofiles  tend  to  pass  through  the  origin.  This  characteristic  appears 
to  be  independent  oheok  on  the  diffusion  analysis. 

Comparisons  of  smoke  diffusion  in  stable  and  unstable  atmospheres  at  Datelan, 
Arizona,  are  aiade  in  Figure  11.  Obseirve  that  the  increased  rate  of  vertical  saoka 

brief  consideration  of  ooBvee^ive  heat  transfer  la  the  ^senee  of  ataaspherie 
thann^  is  given  in  tho  diseassion  section. 
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Figiire  9«  Kbmentura  Eddy  Diffusivity  as  a  Function  of  Height, 
for  Several  Wind  Velocities  at  Riverside 
(Janaary,  1950) 
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Thtrmal  Eddy  Diffusivity,  ft^fhr 


Fisarc  10.  Vertical  ^eraal  Eddy  Diffusivity  Profiles 
M  a  Function  of  Tims  at  Oatslant  Arisons 


(Karch,  191^6) 
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diffusion  Is  In  agreement  with  the  increased  eddy  diffusivlty  at  that  time. 

Figure  12  reveals  day  and  night  thermal  eddy  diffusivlty  profiles  in  the 
lower  three  hundred  feet  for  Manor,  Texas.  The  diurnal  stability  characteristics 
observed  in  the  Arisons  analysis  are  also  noted  in  the  Texas  data. 

C.  Eddy  Diffusion  Analyses  Found  in  the  Literature 

Some  of  the  edd*  diffusion  analyses  found  in  the  literature  are  shown  in 
Figures  13  and  14«  Heat  and  momentum  transfer  diffusivlty  profiles  in  the  lower 
twenty  feet  are  presented  in  Figure  13.  Curves  B  and  C  were  deduced  from 
Sheppard's  shear  stress  and  velocity  data.  Sheppard's  diffusivlty  data  are 
believed  to  fall  above  the  Riverside  data  because  of  the  more  turbulent  state 
of  the  atmosphere.  Again  note  the  important  influence  that  wind  velocity  has 
on  the  magnitude  of  the  eddy  diffusivities. 

A  critical  examination  of  the  diffusion  data  given  in  Figures  13  and  14 
has  not  been  made  because  some  of  the  pertinent  parameters  which  are  necessary 
to  make  such  an  examination  were  not  reported  in  the  literature.  Also,  some  of 
the  methods  of  analysis  were  not  satisfactory. 


Figure  12,  Vertical  Thermal  Eddy  Diffusirity  Profiles 
as  a  Function  of  Time  at  Manor,  Texas 
(August,  194d) 
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Flgurt  13 •  Some  Eddy  Diffusivity  Profiles  in  the  Lower  20  Feet 
of  the  Atmosphere  Whleli  Have  Been  Reported  in  the  Literature 
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Figure  14*  Some  Eddy  Diffusivity  Profiles  in  the  Lower  300  Feet 
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Dzscossxoir 

Inrestlgators  In  the  field  of  mlorometeorology  hare  questioned  the  use  of  the 
approzlnate  eddy  diffusion  equations  (idileh  relate  transfer  rates  to  eddy  dlffuslr 
Itles  and  potential  gradients)  under  extremely  unstable  atmospheric  conditions. 
Consider  an  unstable  atmosphere  In  vhloh  large  thernal  eonreetlon  pattenui  are  In 
motion.  It  seems  that  rertlcal  conrectlTe  heat  transfer  should  be  expressed  as 
the  product  of  a  rertloal  conrectlTe  mass  flov  rate,  a  specific  heat,  and  a  dif¬ 
ference  In  m^  tenqparature  of  rising  and  falling  air  masses  at  a  glren  eleratlon. 
Although  some  studies  of  eonreetlon  patterns  ..xre  been  presented  In  the  llteratnre 
(Beferenee  3^),  atteaqpta  to  solre  the  hydrodynamlo  and  heat  transfer  equations  ■ 

'  •  f. 

alaultaneously  for  the  eonreetlon  cell  systea  should  be  made.  Boverer);  If  the, 
rertlcal  air  Velocities  are  not  described  by  a  regular  eonreetlon  pattern  but 
consist  of  random  fluctuations  as  In  the  case  of  stable,  neutral,  or  sllflfhtly  un¬ 
stable  flov,  rertloal  heat,  mass,  and  momentum  transfer  rates  can  be  expressed  by 
the  eddy  diffusion  equations  (l)  and  (2).  These  relations  hare' satis faotorlly 

•  *  e  ••  • 

e 

been  used  In  duct  flov  systems. 

A  slapl«  heat-Bomsntum  transfer  analogy^  for  the  surface  layer  has  been 
developed  vhloh  relates  some  of  the  pertinent  heat  and  momentum  transfer  rarlables 
that  hare  been  dealt  vlth  la  this  paper.  The  surface  layer  is  postulated  to  eon- 
slst  of  a  laminar  sublayer  and  a  'tuzbulent  layer.  The  foUovlng  table  and  figure 
giro  the  heat  and  momentum  transfer  relations  used  la  the  derlratloa. 


^thls  analogy,  vhloh  Is  expresseA  Im  terns  of  dlnsasloalees  groups  or  moiuli. 
Is  similar  to  the  duct  flov  snaXogy  dereloped  by  Boelter,  Nartlnslli,  ami 
Joaasssn  (isfisronoe  9). 


Height 
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Figure  15,  Velocity  and  Temperature  Notations 

It  1.  po»»lbla  to  wl*.  for  tb.  olr  tempwaturo  Inororanto,  tf.  -  T,) 

(T,  -  V.  fro.  tb.  b«t  no«  «lu.Uon.  for  tb.  looinw  and  turbul«.t  lajar.. 
Ib.a.  t«>,«rat»r.  lnor»«.t.  ««.  b.  Kbatltutad  Into  tb.  follovins  datlning 
■  tor  tb.  unit  tbarSat  tonvMtlr.Tiond«ob«oo  (..fiatlnn  t). - 


vihere,  v'  ■  kinematic  viscosity 

g  ■  acceleration  of  gravity 

Upon  the  substitution  of  the  temperature  increments  (T^  ~  ij)  and  (Ij  ~ 
and  the  dimensionless  moduli  given  by  the  equations  (59)  into  equation  (4)«  the 
follovdng  heat-momentum  transfer  analogy  results t 


It  is  sometimes  desirable  to  express  the  shear  stress  in  terms  of  the  drag  co¬ 
efficient.  ThuSf  equation  (60)  may  be  written  as  follows: 


Equation  (60)  or  (6l)  relates  the  micrometeorological  variables  that  describe  the 
heat  and  momentum  transfer  processes  within  the  surface  layer.  These  variables 
are  the  1)  thermal  conductance^  2)  velx>clty  at  the  reference  height,  3)  drag  co¬ 
efficient,  4)  thickness  of  the  laminar  sublayer,  5)  eddy  diffusivity  ratio,  a'  , 

6)  fl\d4  properties,  and  7)  the  stability  parameters  K  and  sg.  The  moduli  Nu,.  , 

Re^  and  the  drag  coefficient  must  all  be  expressed  in  terms  of  the  same  refel^> 
ence  height,  wherever  it  may  be  in  the  turbulent  layer  (x^  must  be  within  the 
surface  layer) • 

E^qierimental  studies  of  all  of  the  micrometeorological  variables  appearing 
in  equation  (6l)  are  being  conducted  at  tiie  University  of  California.  Some  of 
these  studies  have  been  described  in  this  paper;  other  studies  Involve  the  inves¬ 
tigation  of  the  foUotdng,  under  a  range  of  stability  conditions:  1)  the  laminar 
sublayer  thickness  for  simple  flow  systems,  2)  the  parameter,  a'  ,  from  heat  and 
momentum  transfer  measurements,  jand  3)  the  behavior  of  the  stability  and  rouj^iss 
parameters,  K  and  ^  . 
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